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Abstract The confinement mechanism proposed earlier by the author is applied 
to problem of arising the so-called scale Aqcd within the framework of QCD. 
The natural physical assumption consists of that \ I Aqcd ~ < r > where < r > 
is a characteristic size of hadron (radius of confinement). The above confinement 
mechanism allows us to calculate < r > for mesons in terms of quark and gluonic 
degrees of freedom and this permits to conclude that Aqc d should slightly change 
from hadron to hadron. 
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1 Introduction 



As is known (see, e.g., [I][2]f3] ) . within the framework of perturbative quantum 
chromodynamic (QCD) in asymptotical region of large momenta there appears 
some constant Aqcd with dimension of mass, Aqcd being not connected with 
| arbitrary choice of the renormalization point. This is the so-called scale Aqcd ( m 

what follows we denote it by just A) and it cannot be determined from the theory 
in asymptotical region of large momenta and should arise in some way in the 
region of small momenta. The merits of case is contained in the relation between 
the effective strong coupling constant a s and the momentum transfer \J Q 2 

o: 10 _ 

a s =a s (Q 2 )= (33 _ 2?V )m(Q2M 2 ) W 
while rif is number of quark flavours and (1) holds true at the momentum 
transfer \/Q 2 — > oo. As is clear from (1), if Q 2 — > A 2 then a s — > oo and traditional 
physical interpretation of it consists in that at distances of order r ~ 1/A (r is dis- 
tance between quarks) interaction of quarks becomes very strong and provides the 
quark confinement so at distances greater than 1/A quarks are unobservable (see, 
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e.g., [1,2 ). It should be noted, however, that in fact the numerical value of A de- 
pends on the situation where the relation (1) is used [2]. Under the circumstances, 
appearance of the scale A should be closely connected with the confinement mech- 
anism and the latter should give an evaluation prescription of the mean radius of 
confinement < r > which is in essence just a characteristic size of one or another 
hadron. 

It is known, however, at present no generally accepted quark confinement mech- 
anism exists that would be capable to calculate a number of nonperturbative pa- 
rameters characterizing mesons (masses, radii, decay constants and so on) appeal- 
ing directly to quark and gluon degrees of freedom related to QCD-Lagrangian. 
At best there are a few scenarios (directly not connected to QCD-Lagrangian) 
of confinement that restrict themselves mainly to qualitative considerations with 
small possiblities of concrete calculation. In view of it in [5j[6l[I] a confinement 
mechanism has been proposed which was based on the unique family of compat- 
ible nonperturbative solutions for the Dirac- Yang-Mills system directly derived 
from QCD-Lagrangian. The word unique should be understood in the strict math- 
ematical sense. Let us write down arbitrary SU(3)-Yang-Mills field in the form 
A = A^dx^ = A^Xadx* 1 (A a are the known Gell-Mann matrices, /i = t, r, $, tp, 
a = 1,...,8 and we use the ordinary set of local spherical coordinates r,d,tp for 
spatial part of the flat Minkowski spacetime) . 

In fact in 5,6,7| the following theorem was proved: 

The unique exact spherically symmetric (nonperturbative) solutions (depending only 
on r and r^ 1 ) of SU(3)-Yang-Mills equations in Minkowski spacetime consist of the 
family of form 

V3 r V3 r 

A 3t ^-j=Al = ^±^-(A 1 + A 2 ), 

Ai v =A% + -j=A% = b ir + Bi , A 2v = -A% + -j=A% = b 2 r + B 2 , 



A 3v = — J =A% = -(b 1 +b 2 )r-(B 1 +B 2 ) (2) 



_2_ 

with the real constants a,j , Aj , bj , Bj parametrizing the family. Besides in [5J 
[6lE] it was shown that the above unique confining solutions (2) satisfy the so- 
called Wilson confinement criterion 8,9J. Up to now nobody contested this result 
so if we want to describe interaction between quarks by spherically symmetric 
SU(3)-fields then they can be only the ones from the above theorem. On the other 
hand, the desirability of spherically symmetric (colour) interaction between quarks 
at all distances naturally follows from analysing the pp-collisions (see, e.g., [10] ) 
where one observes a Coulomb-like potential in events which can be identified with 
scattering quarks on each other, i.e., actually at small distances one observes the 
Coulomb- like part of solution (2). Under this situation, a natural assumption will 
be that the quark interaction remains spherically symmetric at large distances 
too but then, if trying to extend the Coulomb-like part to large distances in a 
spherically symmetric way, we shall inevitably come to the solution (2) in virtue 
of the above theorem. 
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The aim of the present paper is to some extent to discuss a possible connection 
between the above confinement mechanism and the appearance of the scale A in 
QCD. Section 2 contains a description of the confinement mechanism in question. 
Section 3 is devoted to how the mean radius of confinement could be computed 
within the framework under discussion while Section 4 is devoted to discussion and 
concluding remarks. Finally, Appendix A supplements the present section with a 
proof of the above uniqueness theorem in the case of SU(3)-Yang-Mills equations. 



2 Quark confinement mechanism 

The applications of the family (2) to the description of both the heavy quarkonia 
spectra |lllll2|[T3lll4j and a number of properties of pions, kaons, 77- and ?/-mesons 
15, 16, 17, 18, 19,20 showed that the confinement mechanism is qualitatively the 
same for both light mesons and heavy quarkonia. At this moment it can be decribed 
in the following way. 

The next main physical reasons underlie linear confinement in the mechanism 
under discussion. The first one is that gluon exchange between quarks is realized 
with the propagator different from the photon-like one, and existence and form of 
such a propagator is a direct consequence of the unique confining nonperturbative 
solutions of the Yang-Mills equations [5l|6j[7]. The second reason is that, owing to 
the structure of the mentioned propagator, quarks mainly emit and interchange 
the soft gluons so the gluon condensate (a classical gluon field) between quarks 
basically consists of soft gluons (for more details see [SlIMI]) but, because of the 
fact that any gluon also emits gluons (still softer), the corresponding gluon con- 
centrations rapidly become huge and form a linear confining magnetic colour field 
of enormous strengths, which leads to confinement of quarks. This is by virtue of 
the fact that just the magnetic part of the mentioned propagator is responsible 
for a larger portion of gluon concentrations at large distances since the magnetic 
part has stronger infrared singularities than the electric one. In the circumstances 
physically nonlinearity of the Yang-Mills equations effectively vanishes so the lat- 
ter possess the unique nonperturbative confining solutions of the Abelian-like form 
(2) (with the values in Cartan subalgebra of SU(3)-Lie algebra) which describe the 
gluon condensate under consideration. Moreover, since the overwhelming majority 
of gluons is soft they cannot leave the hadron (meson) until some gluons obtain 
additional energy (due to an external reason) to rush out. So we also deal with 
the confinement of gluons. 

As has been repeatedly explained in p^ll 1 2 U 1 3 U 1 4iri5lfT6irTTl [T8 lfT9l[20] . parame- 
ters Ai t 2 of solution (2) are inessential for physics in question and we can consider 
Ax = A<2 = 0. Obviously we have Ylj=i >4jt = X!j=i A/v> = which reflects the 
fact that for any matrix T from SU(3)-Lie algebra we have TrT = 0. Also, as 
has been repeatedly discussed by us earlier (see, e. g., [11,12,13,14,15,16,17,18, 
19,20]), from the above form it is clear that the solution (2) is a configuration de- 
scribing the electric Coulomb-like colour field (components A^' 8 ) and the magnetic 
colour field linear in r (components A^ s ) and we wrote down the solution (2) in 
the combinations that are just needed further to insert into the Dirac equation 
(directly derived from QCD-Lagrangian) giving the interaction energy between 
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two quarks in a meson and which will look as follows [after inserting solution (2)] 

(3) 



\dtPa J 



E1> = 



E x 
E 2 
E 3 




E 2 *2 

\ E 3 V 3 J 



where Hamiltonian Ej is 



E H 



■In . 2l ( a . 1 1 2 \ 
fJ.0 - 17 Or - i-y - I d$ + -7 7 I 



«7 



1 



r sin i? 



<9^ + — sini97 1 7 3 



+ — cosi?7 2 7 3 



A 1 3 
-57 ( 7 Ajt + 7 



1 



r sin ■& 



(4) 



with the gauge coupling constant g, /iq is a mass parameter and one should consider 
Ho = m qi m q2 1 (m qi + m q2 ) to be the reduced mass composed of the current masses 
m qi .q 2 of quarks forming a meson (quarkonium) while & = (S^i, #2, ^3) with the 
four-dimensional Dirac spinors 9j representing the jth colour component of the 
meson, so 9 may describe relative motion of two quarks in mesons. 

Additional considerations show that the unique modulo square integrable (non- 
perturbative) solutions of the Dirac equation (3) in the field (2) (i.e. relativistic 
bound states) are (with Pauli matrix a\, for more details see [MM?]) 



9* 



1,2,3 



(5) 



with the 2D eigenspinor <Pj 



'-31 



of the euclidean Dirac operator Vq on the 



unit sphere S 2 , while the coordinate r stands for the distance between quarks. The 
explicit form of <Pj is not needed here and can be found in [1U IT2, 13, 14, 15 16, 
I171I181[T91|20| . One can only remark that spinors &j form an orthonormal basis in 
h\ (8 2 ) (in what follows we denote L%{F) the set of the modulo square integrable 
complex functions on any manifold F furnished with an integration measure, then 
L^{F) will be the n-fold direct product of L2(F) endowed with the obvious scalar 
product). 

We can call the quantity uij relative energy of jth colour component of meson 
(while tpj is wave function of a stationary state for jth colour component) . Under 
this situation, if a meson is composed of quarks with different flavours then 
the energy spectrum of the meson will be given by e = m qi + m q2 + uj with the 
current quark masses m qk ( rest energies) of the corresponding quarks and we 
should put the interaction energy cj = ljj for any j in virtue of Dirac equaion (3). 
On the other hand, ljj is given by (for more details see Refs. [MM?]) 



A j g 2 a j b j ± (rij + aj)yj(n? + 2njctj + A?)/j% + g 2 b 2 (n 2 + 2njQj) 



(6) 
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where a 3 = -(ai + a 2 ), 63 = -(61 + 62), B 3 = -(B± + B 2 ), Aj = \j - gBj, 
ctj = \J A 2 - — g 2 a 2 , nj = 0,1,2,..., while \j = ±(lj + 1) are the eigenvalues of 

euclidean Dirac operator X>o ° n unit sphere with lj = 0, 1, 2, .... 

In line with the above we should have the interaction energy oj = u>x = w 2 = ujs 
in energy spectrum e = m qi + m q2 +uj for any meson (quarkonium) and this at once 
imposes two conditions on parameters ay, bj, Bj when choosing some experimental 
value for e at the given current quark masses m qi , m q2 . 

For example, for reference we shall give the radial parts of (5) at nj = (the 
ground state) that are 

Fji = Cj'','" ' ' ' (l Pi = 9bj + P 3 , 

F j2 = U ■,(),,■■' e~^ r (l + &\ , Qj = Mo - Uj (7) 

with [3j = <J [Iq — oj 2 + g 2 b 2 while Cj is determined from the normalization con- 
dition {\Fji\ 2 + \Fj 2 \ 2 )dr = |. Consequently, we shall gain that &j G L^M 3 ) 
at any t € K and, as a result, the solutions of (6) may describe ground state of 
mesons or quarkonia. The same holds true for wave functions (5) in general form 
so that (5) may describe relativistic bound states of mesons (quarkonia) with the 
energy (mass) spectrum e = m qi + m q2 + to. 



3 Meson spectroscopy and the scale Aqcd 

We would take a sin upon our soul if we did not try using the above unique 
family of compatible nonperturbative solutions for the Dirac- Yang-Mills system 
directly derived from QCD-Lagrangian for description of confinement since such 
an approach to confinement has good mathematical and physical grounds and, as 
a result, the approach is itself unique, nonperturbative and relativistic from the 
outset and it appeals immediately to quark and gluon degrees of freedom as should 
be from QCD-point of view. 

Having the off-the-shelf wave functions for mesons of the form (5) we can 
directly calculate < r > in accordance with the standard quantum mechanics rules 



as < r >= J J r 2 ^^d 3 x = , / ^ J r 2 >P^'Pjd 3 x and the result will, e.g., for the case 
of ground state (7), be equal to (for more details see |141ll5[|16U17[[T8lll91l20] ) 



< r >= 



and it is just the radius of meson (quarkonium) determined by the wave func- 
tions of (7) (at rij = = lj) with respect to strong interaction, i.e., radius of 
confinement. 

As an illustration of (6) and (8), Tables I and II (taken from [19]) contain 
the relevant numerical results for charged pions and kaons where it was accepted 
m u = 2.25 MeV, m d = 5 MeV, m s = 107.5 MeV when computing. 
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Table 1 Gauge coupling constant, reduced mass no and parameters of the confining SU(3)- 
gluonic field for charged pions and kaons 



Particle 


9 


MO (MeV) 


0,1 


«2 


bi (GeV) 


b 2 (GeV) 


Bi 


B 2 


k ± — ud, ud 


6.09131 


1.55172 


0.0473002 


0.0118497 


0.178915 


-0.119290 


-0.230 


0.230 


K ± — us, us 


5.30121 


2.20387 


0.167182 


-0.0557501 


0.120150 


0.131046 


-0.900 


0.290 



Table 2 Theoretical and experimental masses and radii for charged pions and kaons 



Particle 


Thcorct. it (McV) 


Experim. it (MeV) 


Thcorct. < t > (fm) 


Experim. < r > (fm) 


7T ± — ud, ud 


a — rn u + m d + Uj (0, 0, 1) = 139.570 


139.56995 


0.673837 


0.672 


K ± — US, US 


a = m u + m s + uij(0, 0, 1) = 493.677 


493.677 


0.544342 


0.560 



Then, as was said early in the paper, we should put l/A ~ < r > with < r > 
of (8) and it is clear that A is a function of the parameters of the confining SU(3)- 
gluonic field (2) and the current quark masses, i.e., A is expressed through gluonic 
and quark degrees of freedom, as should be according to the first principles of 
QCD. Also it is clear that because of < r > slightly changes from hadron to hadron 
(see Table II) then A should do it as well. It should be noted that crucial role in 
generating < r > or A belongs to the magnetic colour field of (2) (parameters 

bj). Indeed, as follows from (8) at \bj\ — > oo we have < r >~ y ^ TaWTF , SO ' n 

the strong magnetic colour field when |6j| — > oo, < r >— > 0, while the meson wave 
functions of (5) and (7) behave as \Pj ~ e _9 ' &J ' r , i.e., just the magnetic colour field 
of (2) provides two quarks with confinement. Thus, our confinement mechanism 
gives an explanation of appearing the scale A in QCD. 



4 Concluding remarks 

The results of present paper as well as the ones of [TT1|12II13U14U15U16H17I[18H191 
[2"U] allow one to speak about that the confinement mechanism elaborated in [31 
E1E1 gives new possibilities for considering many old problems of hadronic (me- 
son) physics (such as nonperturbative computation of decay constants, masses 
and radii of mesons [15,16,17,18,19,20 , chiral symmetry breaking 19,20 and so 
forth) from the first principles of QCD immediately appealing to the quark and 
gluonic degrees of freedom. This is possible because the given mechanism is based 
on the unique family of compatible nonperturbative solutions for the Dirac- Yang- 
Mills system directly derived from QCD-Lagrangian and, as a result, the approach 
is itself nonperturbative, relativistic from the outset, admits self-consistent non- 
relativistic limit and may be employed for any meson (quarkonium) . Under the 
circumstances the words quark and gluonic degrees of freedom make exact sense: 
gluons come forward in the form of bosonic condensate described by parameters 
aj, bj, Bj from the unique exact solution (2) of the Yang-Mills equations while 
quarks are represented by their current masses m q . 
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Finally, one should say that the unique confining solutions similar to (2) ex- 
ist for all semisimple and non-semisimple compact Lie groups, in particular, for 
SU(iV) with N > 2 and U(JV) with TV > 1 [5;. Explicit form of solutions, e.g., 
for SU(iV) with N = 2, 4 can be found in [7] but it should be emphasized that 
components linear in r always represent the magnetic (colour) field in all the men- 
tioned solutions. Especially the case U(l)-group is interesting which corresponds 
to usual electrodynamics. Under this situation, as was pointed out in [61E] there is 
an interesting possibility of indirect experimental verification of the confinement 
mechanism under discussion. Indeed the confining solutions of Maxwell equations 
for classical electrodynamics point out the confinement phase could be in electro- 
dynamics as well. Though there exist no elementary charged particles generating a 
constant magnetic field linear in r, the distance from particle, after all, if it could 
generate this elecromagnetic field configuration in laboratory then one might study 
motion of the charged particles in that field. The confining properties of the men- 
tioned field should be displayed at classical level too but the exact behaviour of 
particles in this field requires certain analysis of the corresponding classical equa- 
tions of motion. Such a program has been recently realized in [21]. Motion of a 
charged (classical) particle was studied in the field representing magnetic part of 
the mentioned solution of Maxwell equations and it was shown that one deals 
with the full classical confinement of the charged particle in such a field: under 
any initial conditions the particle motion is accomplished within a finite region of 
space so that the particle trajectory is near magnetic field lines while the latter 
are compact manifolds (circles). Those results might be useful in thermonuclear 
plasma physics (for more details see |21|). 



Appendix A 

The facts adduced here have been obained in Refs. [6l[7] and we concisely give 
them only for completeness of discussion. 

To specify the question, let us note that in general the Yang-Mills equations 
on a manifold M can be written as 

d* F = g(*F A A- A A *F) , (A.l) 

where a gluonic field A = A^dx^ = A^Xadx* 1 [\ a are the known Gell-Mann matri- 
ces, fj, = t, r, tp (in the case of spherical coordinates), a = 1, 8], the curvature 
matrix (field strentgh) F = dA + gAAA = F^ v \ a dx^ Ad/ with exterior differential 
d and the Cartan's (exterior) product A, while * means the Hodge star operator 
conforming to a metric on manifold under consideration, g is a gauge coupling 
constant. 

The most important case of M is Minkowski spacetime and we are interested 
in the confining solutions A of the SU(3)-Yang-Mills equations. The confining 
solutions were defined in Ref. [5] as the spherically symmetric solutions of the 
Yang-Mills equations (1) containing only the components of the SU(3)-field which 
are Coulomb- like or linear in r. Additionally we impose the Lorentz condition on 
the sought solutions. The latter condition is necessary for quantizing the gauge 
fields consistently within the framework of perturbation theory (see, e. g. Ref. [22]), 
so we should impose the given condition that can be written in the form div(yl) = 0, 
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where the divergence of the Lie algebra valued 1-form A = A^dx^ = A^Xadx^ is 
defined by the relation (see, e. g., Ref. [23]) 

divA = *{d*A) = -^dniVSg^ A u ) . (A.2) 

It should be emphasized that, from the physical point of view, the Lorentz con- 
dition reflects the fact of transversality for gluons that arise as quanta of SU(3)- 
Yang-Mills field when quantizing the latter (see, e. g., Ref. [22]). 

We shall use the Hodge star operator action on the basis differential 2-forms 
on Minkowski spacetime with local coordinates t, r, p in the form 

*(dt A dr) = —r 2 sintfdtf A dip , *(dt A d#) = sini9dr A dp , 

*(dt A dp) = ^-^dr A dfi , *(dr A dfi) = sin#Gte A dp , 

smv 

*(dr A dp) = — dt A dtf , *(d$ A dtp) = dt A dr , (A.3) 

smi9 Hsmw 

so that on 2-forms * 2 = —1. More details about the Hodge star operator can be 
found in [23] . 

The most general ansatz for a spherically symmetric solution is A = At{r)dt + 
A r (r)dr + A${r)d-d + A v (r)dp. But then the Lorentz condition (A.2) for the given 
ansatz gives rise to 

sin $dr{r 2 A r ) + a^(sin^) = 0, 

which yields A r = ^ — c °%^ J A$(r)dr with a constant matrix C. But the confining 
solutions should be spherically symmetric and contain only the components which 
are Coulomb-like or linear in r, so one should put C = A$(r) = 0. Consequently, 
the ansatz A = At(r)dt + A v (r)dp is the most general spherically symmetric one. 

For the latter ansatz we have F = dA + gA f\A = —d r Atdt Adr + d r A v dr Adp> + 
g[At, A v ]dt A dp, where [■,•] signifies matrix commutator. 

Then, according to (A.3), we obtain 

*F = (r 2 sin^dr-Atdfi A dp — d r A v dt A dd- 

smw 



sini? 

which entails 



[A t ,A v ]dr Ad#, (A.4) 



while 



d* F = s\w&dr{r 2 d r At) dr A d$ A dp + -J—rd 2 A v dt A dr A d-d , (A. 5) 

sm.it 



*F A A - A A *F = 

r 2 sm-£>\d r A t , A t ] J-^ldrA^, A J \ dt Adfi Adp 

sin tf ' 



smi 



[LAtjAp], A t ] dt A dr A di3 



+ [[A t ,A<p],A v ] dr A dtf A dp I . (A.6) 
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Under the circumstances the Yang-Mills equations (A.l) are tantamount to the 
conditions 

dr (r 2 d r A t ) = — 3[[A t ,A v },A v ], (A.7) 
sin v 

dtA v = -g 2 [[At,A v ],A t ], (A.8) 
r 2 sm^[d r At,A t ] - -J— [d r A v ,A v ] = 0. (A.9) 

The key equation is (A.7) because the matrices At,A v depend on merely r and 
(A.7) can be satisfied only if the matrices At = AfX a and A v = A^\ a belong to the 
so-called Cartan subalgebra of the SU(3)-Lie algebra. Let us remind that, by defi- 
nition, a Cartan subalgebra is a maximal abelian subalgebra in the corresponding 
Lie algebra, i. e., the commutator for any two matrices of the Cartan subalgebra 
is equal to zero (see, e.g., Ref. [23]). For SU(3)-Lie algebra the conforming Cartan 
subalgebra is generated by the Gell-Mann matrices A3, As which are 



(A.10) 



Under the situation we should have At = At A3 + Af \g and A v = A^X^ + A^Xs, 
then [At, A v ] = and we obtain 

d r {r 2 d r At) = 0,drA v = 0, (A.ll) 

while (A.9) is identically satisfied and (A.ll) gives rise to the solution (2) with 
real constants aj, Aj, bj, Bj parametrizing the solution which proves the uniqueness 
theorem of Section 1 for the SU(3) Yang-Mills equations. 
The more explicit form of (2) is 

A 3 t = [{a 2 -a 1 )/r + A 1 ~A 2 }/2, 

Al = [Ai +A 2 - (ai + a 2 )/r]v / 3/2 , 

A% = [(bi-b 2 )r + B 1 -B 2 ]/2, 
A% = [(bi + 6 2 )r + Bi+B 2 ]v / 3/2. (A.12) 

Clearly, the obtained results may be extended over all SU(iV)-groups with N > 
2 and even over all semisimple compact Lie groups since for them the corresponding 
Lie algebras possess just the only Cartan subalgebra. Also we can talk about the 
compact non-semisimple groups, for example, U(iV). In the latter case additionally 
to Cartan subalgebra we have centrum consisting from the matrices of the form 
a/jv (In is the unit matrix iV x N) with arbitrary constant a. 

The most relevant physical cases are of course U(l)- and SU(3)-ones (QED 
and QCD). In particular, the U(l)-case allows us to build the classical model of 
confinement (see Ref. [21]). 

At last, it should also be noted that the nontrivial confining solutions obtained 
exist at any gauge coupling constant g, i. e. they are essentially nonperturbative 
ones. 



10 



Yu. P. Goncharov 



References 



1. K. Huang, Quarks, Leptons and Gauge Fields (World-Scientific, Singapore, 1982). 

2. F. J. Yndurain, Quantum Chromodynamics. An Introduction to the Theory of Quarks and 
Gluons (Springer- Verlag, Berlin, 1983). 

3. D. W. Duke, R. G. Roberts, Phys. Rep. 120, 275 (1985). 

4. Particle Data Group (K. Nakamura, et al), J. Phys. G: Nucl. Part. Phys. 37, 075021 (2010). 

5. Yu. P. Goncharov, Mod. Phys. Lett. A 16, 557 (2001). 

6. Yu. P. Goncharov, Phys. Lett. B 617, 67 (2005). 

7. Yu. P. Goncharov, in: P. V. Kreitler (Ed.), New Developments in Black Hole Research, pp. 
67-121. Nova Science Publishers, New York (2006), Chap.3, [hep-th/0512099| 

8. K. Wilson, Phys. Rev. D 10, 2445 (1974). 

9. M. Bander, Phys. Rep. 75, 205 (1981). 

10. D. H. Perkins, Introduction to High Energy Physics (Cambridge University Press, Cam- 
bridge, 2000). 

11. Yu. P. Goncharov, Europhys. Lett. 62, 684 (2003). 

12. Yu. P. Goncharov, E. A. Choban, Mod. Phys. Lett. A 18, 1661 (2003). 

13. Yu. P. Goncharov, A. A. Bytsenko, Phys. Lett. B 602, 86 (2004). 

14. Yu. P. Goncharov, Nucl. Phys. A 808, 73 (2008). 

15. Yu. P. Goncharov, Phys. Lett. B 641, 237 (2006). 

16. Yu. P. Goncharov, Phys. Lett. B 652, 310 (2007). 

17. Yu. P. Goncharov, Mod. Phys. Lett. A 22, 2273 (2007). 

18. Yu. P. Goncharov, J. Phys. G: Nucl. Part. Phys. 35, 095006 (2008). 

19. Yu. P. Goncharov, Nucl. Phys. A 812, 99 (2008). 

20. Yu. P. Goncharov, Eur. Phys. J. A 46, 139 (2010). 

21. Yu. P. Goncharov, N. E. Firsova, Int. J. Theor. Phys. 49, 1155 (2010). 

22. L. H. Ryder, Quantum Field Theory (Cambridge Univ. Press, Cambridge, 1985). 

23. A. L. Besse, Einstein Manifolds (Springer- Verlag, Berlin, 1987). 

24. A. O. Barut, R. Raczka, Theory of Group Representations and Applications (Polish Sci- 
ence, Warszawa, 1977). 



